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Abstract 

Hexagonal polyominoes are polyominoes on the honeycomb lattice. We enumerate the 
symmetry classes of convex hexagonal polyominoes. Here convexity is to be understood as 
convexity along the three main column directions. We deduce the generating series of free (i.e. 
up to reflection and rotation) and of asymmetric convex hexagonal polyominoes, according 
to area and half-perimeter. We give explicit formulas or implicit functional equations for the 
generating series, which are convenient for computer algebra. 

1 Introduction 

An hexagonal polyomino is a finite connected set of basic cells of the honeycomb lattice in the 
plane. Note that the hexagons of our lattice have two sides parallel to the horizontal axis. See 
Figure ^ Unless otherwise stated, all the polyominoes considered here are hexagonal. The 
area of a polyomino is the number of the cells composing it. Its perimeter is the number of line 
segments on its boundary. We say that a polyomino is convex along a direction if the intersection 
with any line parallel to this direction and passing through the center of a cell is connected. 
The convexity directions are characterized by the angle a (0 < a < vr) which they form with 
the positive horizontal axis. 

Various convexity concepts have been introduced in the literature for hexagonal polyominoes, 
depending on the required convexity directions. Following the nomenclature of Denise, Diir, 
and Hassani [2j, we mention the EG-convex polyominoes, where a = and '/r/2, studied by 
Guttmann and Enting [5] and by Lin and Chang |^, the C^-convex polyominoes, where a = 7^/2, 
enumerated according to many parameters by Lin and Wu ^Hl and by Feretic and Svrtan [1], 
the strongly convex polyominoes, where a = 0, tt/3 and 27r/3, introduced by Hassani |H] and 
studied in [£] and |2j, and finally the C- or -convex polyominoes, where a = vr/G, 7r/2 and 
57r/6, introduced and enumerated according to perimeter in jB] and In particular, Hassani 
gives explicitly the algebraic generating function for C-convex polyominoes according to half 
perimeter. 

It is this last class that interest us here, and that we call simply convex polyminoes. See 
Figure n for an example. This concept is a natural extension of (row and column) convexity on 
the square lattice. 

*With the partial support of CNRS (France), NSERC (Canada) and FCAR (Quebec). This is the full version 
of a paper presented at the FPSAC Conference in Vancouver, Canada, June 28 - July 2, 2004 (see 0). 
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Figure 1: A convex polyomino and its convexity directions 



These polyominoes are traditionally taken up to translation. However it is natural to con- 
sider them also up to rotation and reflection, as objects living freely in space. Following Voge, 
Guttmann and Jensen , we call these equivalence classes free polyominoes. In organic chem- 
istry, free polyminoes represent benzenoid hydrocarbons. See where these molecules (with- 
out the convexity property) are enumerated by an exhaustive generation approach. 

Our objective is to enumerate free convex polyominoes, according to area and half perimeter. 
Following the approach of Leroux, Rassart and Robitaille |H] for the square lattice, we consider 
them as orbits of the dihedral group Vq (the group of isometrics of the regular hexagon) , acting 
on convex polyominoes. The Cauchy-Frobenius Formula (alias Burnside's Lemma) can be used 
to count the orbits and we are thus lead to enumerate the symmetry classes Fix(/i) of convex 
polyominoes, for each element h of the group Pg; 

It is also possible to enumerate convex polyominoes which are asymmetric or which have 
exactly the symmetries of a given subgroup H Vq, using Mobius inversion in the lattice 
of subgroups of Pe- For this purpose, we also enumerate the symmetry classes F^h for each 
subgroup H oi Vq. 

For any class JF of convex polyominoes, we denote its generating series by T{x,q,u,v,t), 
where the variable x marks the number of columns, q marks the area, u marks the size of the 
first column (on the left), v, the size of the last column, and t the half perimeter; for example, 
the polyomino of Figure ^ has weight x^^q^^u^v^t^^ . It is possible that some variables do not 
appear in some generating series. The generating series will be given by explicit formulas or 
implicit functional equations, which are convenient for computer algebra. 

2 Preliminaries 

2.1 Particular classes of convex polyominoes 

Some familiar classes of convex polyominoes of the square lattice are naturally found on the hon- 
eycomb lattice and are useful. It is the case notably of partition and staircase (or parallelogram) 
polyominoes and of a variant of stack polyominoes. 
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a) b) c) 

Figure 2: Partition and staircase polyominoes 



2.1.1 Partition polyominoes 

Figure[2K represents the partition (4, 2, 2) contained in a rectangle of size 5 x 4 in the honeycomb 
lattice. Figure ISb represents the distinct part partition (5, 4, 2, 1), with parts bounded by 6. We 
denote by Dm{u, q) the generating polynomial of distinct part partitions with parts bounded by 
m. Here the variable u marks the number of parts. We have 

Dm{u,q) = {l + uq){l+uq^)---{l + uq"') and Do{u,q) = l. (1) 



2.1.2 Staircase polyominoes 

Figure |2t represents a staircase polyomino from the square lattice (see for example ^ or [Jj) 
redrawn on the honeycomb lattice. Observe that the half perimeter is equal to 2p — 1 where p 
is the half perimeter on the square lattice. We know that these polyominoes are enumerated 
according to half perimeter by the Catalan numbers and according to area by the sequence 
M1175 of m] (A006958 of JH]) whose generating series is a quotient of two g-Bessel functions. 

We denote by Pa, the set of staircase polyominoes (Pa for Parallelogram) on the honeycomb 
lattice and by Pa{x,q,u,v,t), their generating series. An analysis of the situation where a 
column is added on the right, following the method of M. Bousquet-Melou (compare with 
Lemma 3.1), gives, for Fa{v) = Fa{x, q,u,v,t), 

Pa(^;) = + ^ (Pa(l) - Pa{vq)) (2) 

1-quvt^ (1 - qvt^){l - qv) ^ v ; v -i;; \J 

so that 

p . ^ _ Jl{l)+Jl{v)Jo{l)-Jl{l)Mv) 

Pa(^) - J^) ' 

where 

h iqvt^;q)n{qv;qUl-q^+'uvt^) 

and 

''^""^-y-'^ iqvt^;qUqv-q)^- 
Here we have used the familiar notation (a; q)n = (1 — a)(l — aq) • • • (1 — aq^~^. We set 

Pa(x, g, u, f , t) = ^ Faij{x,q,t)u^v^ . (4) 
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2.1.3 Stack polyominoes 

There exists a specific variant of stack polyominoes on the honeycomb lattice. They consist of 
pyramidal stackings of hexagons, viewed sideways for our purposes. A first class (see FigurelSK), 
denoted by T (for French tas), appears in the literature under the name of pyramidal stacking 
of circles] see Their generating series according to area is referenced as number M0687 in 
[111 and A001524 in 





b) ^ c) 

Figure 3: Stacks and symmetric stacks 

Let T(x, n, q) be the generating series of stack polyominoes according to the number of 
columns (the width, marked by x), the size of the first column (the height, marked by u) and 
the area, and let Tn{x,q) = [u"]T(x,m, g) be the generating series of stacks whose first column 
is of size n. Note that the half perimeter is equal to twice the height plus the width so that the 
series T{xt,ut^ ,q) also keeps track of this parameter. 

We have 

T{x,u,q)=Y. Y J; (5) 

and 

r„(x, q)=Y. X-q-H-) ^ [--If] ^ . (6) 

m=l j=0 

The polynomials Tn{x,q) can also be rapidly computed by recurrence using the class T0„ of 
stack polyominoes whose first column is of size n, including empty cells at the two extremities. 
See Figure 03. Indeed, we have 

r„(x,(?) =xg'^TO„_i(x,g). (7) 

with TOo(x,g) = 1, TOi(x,g) = 1 + xq, and, arguing on the existence of empty cells at each 
extremity, 

TO„(x, g) = (2 + xg")TO„_i - T0„_2, (8) 

2.1.4 Symmetric stacks 

Horizontally symmetric stacks (see Figures ISt andOJi), constitute the families TS and TSO. 
Using the same notation as for stacks, we have 

^ ^ x^u'^q^^'^+^y^{l+uq^) 

l^[X, U, q) _ _ ^2^4) . . . (1 _ ^2^2m^ ' 
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Moreover, 

TSn{x,q)=xq''TSOn-i{x,q). (10) 
where TSOo(x,g) = TSO_i(a;,g) = 1, and 

TSOn(x, q) = x(?"TSO„_i(x, q) + TS0„_2(x, q), (11) 

2.2 The dihedral group Vq 

The dihedral group "Dq is defined algebraically by 

Vq = {p,t\p^ = 1, = 1, rpr = p~^). 

Here is realized as the group of isometries of a regular hexagon, with p = r = the (clockwise) 
rotation of 7r/3 radian and r = dss = h, the horizontal reflection. We have 

Vq = {id,r,r^,r^,r^,r^,dai,da2,da3,dsi,ds2,ds3}, 

where ds2 = rr^, dsi = rr^, reflections according to vertex- vertex axes, and das = rr, da2 = rr^, 
and dai = rr^, reflections according to the edge-edge axes. See Figure 0] 

dsi da2=v 













^^^^^^ ' 







ds3=h 



Figure 4: The reflections of Vq 

The dihedral group Vq acts naturaly on (hexagonal) polyominoes, by rotation or reflection. 
For any class J- of polyominoes, with a monomial weigh w corresponding to certain parameters, 
we denote by \ J^\w the total weight (i.e. the generating series) of this class. If ^is invariant under 
the action of Vq, then the set of orbits of this action is denoted by J- /Vq. Burnside's Lemma 
enumerates these orbits in terms of the sets Y\yi[g) of fixed points of each of the elements g 
of Vq, the symmetry classes of J-. We write fix((7) = |Fix(5()|^^,. Clearly we have fix(r) = 
fix(r^), fix(r^) = fix(r^) and, for symmetry reasons, fix(dai) = fix(da2) = fix(da3) and fix(dsi) = 
fix(ds2) = fix(ds3). In the following, we will choose v = da2, the vertical axis, and h = dsa, the 
horizontal axis. We then have 

1-^/^6 U = T^Y^ fix(5() 

= ^(|^U + 2fix(r) + 2fix(r2)+fix(r3)+3fix(v) + 3fix(h)). (12) 
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2.2.1 The lattice of subgroups of Pg 



It is also possible to enumerate the convex polyominoes which are asymmetric or which have 
exactly the symmetries of a given subgroup H of T^q, with the help of Mobius inversion in the 
lattice of subgroups of Pg. This lattice and its Mobius function are well described in Stockmeyer's 
Ph.D. thesis ^Hl) for any dihedral group P„. We follow his nomenclature. Apart from the trivial 
subgroups = {id} and 1 = Vq =< ic,]i >, we have the cyclic subgroups 

Cg =< r >= {l,r,r^r^r^r5}, =< r^ >= {1, r^, r^} and ^2 =< r^ >= {l,r3}, 

=< ds2 >= {l,ds2}, Fi,2 =< dsi >= {l,dsi} and Fi,3 =< h >= {1, h}, 
-ffi^i =< da3 >= {1, das}, ^^i,2 =< v >= {1, v} and i^i^g =< dai >= {1, dai}, 
as well as the two generator subgroups 

F3,i =< r2,ds2 >= {l,r^r^dsl,ds2,h} =< r^h >, 

H-i^i =< r2,da3 >= {1, r^, r^, dai, v, daa} =< v'^,v >, 
and the D2J =< r^, rr^-' >, j = 1, 2, 3, that is 

L»2,i =< r^ds2 >= {I,r3,ds2,dai}, L>2,2 =< r^dsl >= {1, r^ dsi, das}, 

and 1)2,3 =< r^,h >= {l,r^,h,v} =< r^,h > . 
The lattice of subgroups of Pg is represented in Figure |SJ 




(-1) 



and 



0={1) 



Figure 5: The lattice of subgroups of Dq (/u(0, H) in parenthesis) 
For any subgroup H of Pg (H < Pg ), we set 

F>H = \{seJ'\ stab(s) > H}\^ = \{s£j^\h£H^h-s = s}\^ 

F=H = \{seJ^\ stab(s) = H}\^ = \{s£j^\heH^h-s = s}\^ 



(13) 
(14) 
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We have clearly, for any H <T>q 

F>H = E 

H<K<V6 

and, by Mobius inversion, 

H<K<Ve 

In particular, the total weight of asymmetric polyominoes is given by 

F=o= E f^iO,K)F>K (15) 

K<V6 

Note that F>o = \J^\w and that for any cyclic subgroup H =< h >, F^h = fix(/i). For reasons 
of symmetry (or by conjugation), we have -F>_D2 i = -^>D2 2 — F>d2z- -^-^ following, we will 
take D2,3 =< r^,h >. The formula (|15|) then yields 

F=o = |^U-3fix(h)-3fix(v)-fix(r2)-fix(r3)+fix(r)+6F>fl,3+3F>i.3,+3F>H3,i-6F>i,g. (16) 

For any subgroup H of Pgj we sometimes write |Fix(i?)|g^t = F>_ff when the weight is given by 
the area and the half perimeter. 

2.3 Growth phases of convex polyominoes 

Any convex polyomino can be decomposed into blocks according to the growth phases, from left 
to right, of its upper and lower profiles. Figure El gives an example of this decomposition. The 
upper profile is represented by the path from A to B along the upper boundary, and the lower 
profile, by the path from C to D. On the upper profile, we speak of a weak growth if the level 
rises by a half hexagon only with respect to the preceding column, and of a strong growth if the 
level rises by more than a half hexagon. We define analogously a weak or strong decrease. On 
the lower profile, a growth corresponds to a descent and a decrease, to a rise. 




Figure 6: Growth phases of a convex polyomino 



The state in which a column lies is described by an ordered pair (i, j), i, j = 0, 1, 2; the first 
component corresponds to the upper profile and the second, to the lower profile. The state 
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corresponds to a (weak or strong) growth, at the start of the polyomino, the state 1 to a weak 
growth or decrease, in an oscillation phase, and the state 2, to a strong or weak decrease, in 
the last part of the polyomino. To pass from the state to the state 1, there must be a first 
weak decrease, and to pass from the state or 1 to the state 2, there must be a strong decrease. 
Finally, the transitions from the state 1 to the state and from the state 2 to the state 1 or 
are impossible. Now, a block Hij is characterized by a maximal sequence of consecutive columns 
which are in the state 

We can then view a convex polyomino as an assemblage of blocks and we will first enu- 
merate these blocks Hij. In what follows, we give the various generating series of the form 
Hij{x,q,u,v,t). 

2.3.1 The families Hqq and H22 

The polyominoes of the classes Hqq and H22 are easy to enumerate for they are in fact stack 
polyominoes. Here, only one of the two variables u and v is used at a time. We have 



where T{x, u, q) is given by (jSJ. 

2.3.2 The families Fqi, i^io, H12, and H21 

The classes -ffoii -f^iO) H12, and H21 of polyominoes are in bijection with each other by horizontal 
and vertical reflections and are thus equinumerous. Figure [3 shows a polyomino in Hiq. We 
easily find that 



H22{x,q,u,t) = T{xt,ut^,q) and HQQ{x,q,v,t) = T{xt,vt'^,q) 



(17) 



Hio{x,q,u,v,t) 



xquvt^ xt'^{l + qvt) 



Hio{x,q,u,vq,t) 



(18) 






^ x"'q"'uvt'^"'+^{-qvt; q)m-i 
\ {qvt'^;q)m-i{^ - q^uvt"^) 



(19) 




Figure 7: Polyomino of Hiq 
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Formula ()19() can be seen directly on Figure [Tj We also see that 



m— 1 



h>l m>l 



1=0 



'■m-2+j 
j 



j>0 



(20) 



Note that Hoi{x,q,u,v,t) = HiQ{x,q,u,v,t) and that Hi2ix,q,u,v,t) = H2i{x,q,u,v,t) 
Hio{x,q,v,u,t). 



2.3.3 The families H02 and H20 

These two classes are in fact equivalent to staircase polyominoes: 

Ho2ix, q, u, V, t) = Fa{x, q, u, v, t) = H2o{x, q, u, v, t). 



(21) 



2.3.4 The family Hn 



The class Hn contains the convex polyominoes whose upper and lower profiles are both oscillat- 
ing. When we examine the diagonal row of hexagons in the das axis (see Fig. 0} containing the 
first column's lower cell, two subclasses of Hu appear. In the first class, denoted by Hua, this 
diagonal row and those to its right (up to the last column) form a staircase polyomino (rotated 
clockwise by a 7r/3 angle); see Figure |Hh- In the second class, denoted by Hm,, this diagonal 
is the basis of a rectangle of height at least 2; see Figure ISb. In both cases, we find, above 
and below these objects (staircase or rectangle), distinct part partitions which are left and right 
justified, respectively. 





a) b) 
Figure 8: Polyominoes of Hn 

Recall that in the series Paij{x,q,t), defined by @, the variable x marks the number of 
columns of the (unrotated) staircase polyomino. We rather use the link between its width i, 
when rotated, and its half-perimeter p: p = 21 + 1. Hence we have 



Hii{x, q, u, V, t) = Hiia{x, q, u, v, t) + Hubix, q, u, v, t) 



(22) 



with 



and 



Hua{x,q,u,v,t) = X 2uvPaij{l,q,tx'2)Di_i{ut,q)Dj^i{vt,q) (23) 

i>l j>l 

x-q^-uW^+^Dn-ijut, q)Dn-iivt, q) 
Hubix, q, u, v,t) = ^ 1 _ n^^^2 • (24) 
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3 Convex and directed convex polyominoes 

3.1 Convex polyominoes 

We denote by C, the class of all convex polyominoes and by Cij, the subclass of polyominoes 
whose last column is in the state i,j = 0, 1, 2. This determines a partition of C. To enu- 

merate C, we must enumerate each of the classes Cij. We give the generating series Cij{x, q, v, t), 
using the growth phase decomposition of a convex polyomino, following essentially the method 
of Hassani |2] . 

We use the notation dj^Hiiji for the set of convex polyominoes obtained by gluing together 
in all the possible ways a polyomino of Cij with one of Hiiji. We introduce the series Cij,n{x, q, t) 
and Hij^nix, q,v,t) by the coefhcient extractions 

Cij^nix,q,t) = [v"']Cij{x,q,v,t) and Hij^n{x,q,v,t) = [u'^]Hij{x,q,u,v,t). (25) 

For example, we have Cqo = Hqq, Ciq = Cqo ^ Hiq and 

Cio{x, q,v,t) = XI ( X] 72E^'^ooA^^ 1^ ) Hw^x, q, v, t) 

n>l \k=l J 

= J2{j2^^k{xt,qUH^oA^^Q^v^t) (26) 

n>l \fc=l / 

= Coi{x,q,v,t). 

Likewise, Cn = (Cqo + Cio + Coi) ® Hn = Coo ® Hn + Cio ® Hn + Cqi Hn and 

Coo^ Hii{x,q,v,t) = '^-^^Coo,n{^,q,t)Hii.n-i{x,q,v,t), (27) 

n>2 ^ 



1 1 

Cio®Hii (x, q, v,t) = ^ 72;i3TCio,n(2;, g, t)i?ii,„(a;, g, v, i)+X Ti;!^^ Cio,n(a;, q, t)Fii,„-i(a;, g, t). 

(28) 



n>l n>2 



We have also C02 = (Cqo + Cqi) O i?'o2, 

C12 = (Coo + Coi + Cio + Cii + C02) ® H12, 
C22 = (Coo + Coi + Cio + Cii + C02 + C02) (8) -f^22- 

Finally, 

C{x,q,v,t) = {Coo + 2Cio + Cn + 2Co2 + 2Ci2 + C22){x,q,v,t). (29) 



3.2 Directed convex polyominoes 

There is a special class of convex polyominoes which will be particularly useful in the following, 
namely those which are directed according to the North direction with a diagonal basis. See 
Figure ini We call them simply directed convex, and denote this class by D. 

Analogously to convex polyominoes, any polyomino in D can be decomposed into blocks Hi 
according to the growth phases i = 0,1, ou 2, of its upper profile. Their generating series can 
be computed directly by observation. Since the half perimeter can be deduced from the other 
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Figure 9: Directed convex polyomino 



parameters, the variable t does not appear here. For example, a polyomino in Hq is identified 
with a partition and we find 



1>1 k>l 



l+k-2 



A polyomino in Hi can be decomposed into pieces as shown in Figure ITIIk . yielding 

(m)(m+l) r. 



Hi{x,q,u,v) =^v'- ^ u'-'^"^ ^ x'-fl"" ' 2 



l>l m>0 k>m+l 



k~l 
m 



(30) 



(31) 



Likewise, for H2, we find (see Figure ITUb) 





Figure 10: Directed convex polyominos in Hi and in H2 



I I 1^,1 , ^fc ^i+A:+m-l„-!ii^i2^+m 

1>1 \ k>2 m>0 



H2{x, q, n, f ) = ^ f M xq'-u^ + ^ x'^ ^ " ' ^q 2 



m+k-2 
k-2 



(32) 



We denote by Pj, the class of directed convex polyominoes whose last column is in the state 
i, i = 0, 1, 2, and we introduce the notation 



T>i^n{x,q,t) = [v"-]'Di{x,q,v,t) and Hi^n{x,q,v) = [u"-]Hi{x,q,u,v). 



(33) 
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We have 

Vo{x, q, V, t) = -Ho{xt^,q, vt^), (34) 

Vi = Vq® Hi and 

Vi{x,q,v,t) = ^Vo^„i+i{x, q,t)Hi^„i{xt^ , q,vt), (35) 

m>l 

and finally, V2 = (Pq + 'Di) ® H2 and 

V2{x, q,v,t) = "YYl ^{Vo,m+h{x, q, t) + Vi.m+h{x, q, ^))i^2,m(x^^ q, vt). (36) 

m>l h>2 



4 Reflexive symmetry classes 
4.1 Vertical symmetry 

Consider a vertically symmetric (v-symmetric) convex polyomino P. We see that the symmetry 
axis goes through a central column. Denote by K the left fundamental region of P, including 
the central column. See Figure ITTl We have 




Figure 11: Polyomino convex w-symmetric 



K{x,q,v,t) = Coo{x,q,v,t) +2Cio{x,q,v,t) + Cii{x,q,v,t) (37) 
= ^ K^(x,(?,t)i;- 

m>l 

and 

|Fix(v)|,, = Y: Z^^K„,{l,q^t'). (38) 

m>l ^ ^ 

4.2 Horizontal symmetry 

The class S of h-symmetric convex polyominoes is partitioned into three classes: Sa and Sb, 
wether or not we can find an arrowhead polyomino in the oscillating part (see Figures 112b . and 
112b ) and the class Sc, if there does not exist an oscillating part. 

In order to construct a polyomino of the class arrrowhead, denoted by A, we start with a 
triangle of side n, to which a symmetric stack is possibly attached to form the H22 phase; denote 
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Figure 12: /i-symmetric convex polyominoes 



by B, this starting class of polyominoes. From B, we construct A by successively attaching T-- 
shaped bands on the left, as illustrated in Figure [T^ . Let the variable s mark the size of the 
last attached F's upper left part. We have 

B{s, X, q, t) = sxqt^ + ^ s"x"g"("+i)/2^=^nrpgQ^_^(^^^ ^gg^ 

n>2 

and the generating series A{s) = A{s, x, q, t) is characterized by the following functional equation, 
which can be solved by the usual method: 

2 2 3 4 

A{s) = B{s) + ^p^^^(l) - Aisq""). (40) 
1 — sg^ 

We set A{s,x,q,t) = J2m>i^rn{x,t,q)s"^. To complete the polyomino, we must take into 
account the parity of the first oscillating column. The first case, illustrated in Figure 112b .. is 
when this column is of odd size. In the second case, this even column is placed in front of the 
arrrowhead. In conclusion, we obtain 

Sa{x,t,q) = 5](?'^('^+l)t2h+2TS2,+2(^t,g) ^ ^A^{x,t,q) 
h>0 m>h+l ^ 

+ ^gM/^+i)t2/.+3TS2,+i(:rt,g) J2 ['^],2AUx,t,q). (41) 

h>l m>h 

The computations for Sb{x, t, q) and 5c(x, t, q) are simpler. For Sh, there are also two parity 
cases and we find directly 

n-l 

S,{x,t,q) = ^X>(T)i3n^^2fcn^4/=TS0„+2,_3(xt,(?)5]t2h+2^'^(h+l) Jn^l] ^ TS2fc+2/.+2(x, t, g) 
n>l k>l h=0 

n 

+ ^"9^"^')*'" E 9''("+')i''+'TS0„+2fc-3(^t, q) E t^^^^^^^^'^ TS2fc+2h+i(x, t, q) 

n>0 k>l h=0 

(42) 
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and 

Sc{x, t,q) = Y, t^'^TShixt, q)TSOh^s{xt, q). (43) 

h>l 

Finally, 

|Fix(h)|,,j = Saih t, q) + Sb{l,t, q) + 5,(1, t, q). (44) 

5 Rotational symmetry classes 

5.1 Symmetry with respect to the 7r/3 radian rotation r 

The polyominoes which are symmetric with respect to the tt/3 rotation (r-symmetric) are es- 
sentially formed of large hexagons decorated by stack polyominoes of the class TO. We find 

|Fix(r)|<j,j = ^t3{2'^-i)g3M/^-i)+iT0ft_^(t6^g6)^ (45) 

h>l 

the series TO„,(x,g) being defined by (jSJ. 

5.2 Symmetry with respect to the 27r/3 radian rotation 

The situation is more complex here. First we must distinguish the case where the rotation 
center is in the middle of an hexagon from the one where it is on a vertex. This determines two 
subclasses, denoted by V and Q. 

5.2.1 The rotation center is in the middle of an hexagon 

When the rotation center is in the middle of an hexagon, we consider the fundamental region 
formed by the upper third of the r^-symmetric polyomino. The parameters hi and /12 are 
defined as the extent of the fundamental region in the directions da2 = v and dai (or das), 
respectively. There are three subcases according to wether /ii > /i2, ft,2 > hi or hi = /i2, giving 
three subclasses denoted by Vi, V2 and V3. Figure El illustrates the first subcase Vi. 




Figure 13: Fundamental region of an r^-symmetric convex polyomino in Vi 
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In this figure, there is a basis formed of (one third of) a super-hexagon of radius h = h2, on 
top of which are placed a directed convex polyomino D (see the subsection I3.2() and the image 
D* under v of another directed convex polyomino, these two polyominoes sharing a common 
column. Let i be the size of this common column, so that hi = h + L We must consider all 
such legal combinations T) ®T)* and take into account the added area and half perimeter over 
that of the super-hexagon. For i = 0, 1, 2, we set 

A,r,n(g,t) = ■^2^^—[W]Vi^n{x,q,t) = (x, g, ^, t) 

and Vr^niq, t) = (I'o.r.n + 2?i,r,n + T^2,r,n){q, t). We then have 

r,(q^t) = ^t3(2^-l)^3h(h-l)+l 
h>l 

h h 

EEE^ (^o,r,K'?',i')^o,M('?',i') + 22^o,r,/(g',t')2?i,M(9',t')+ 

l>lr=ls=l ^ 

2Po,r,K'^^ ^^)2^2,s,K'7^ i^) + Vi^rAQ^t^)'^i,sAQ^t'^)) ■ (46) 

For reasons of symmetry (a it/3 rotation), we see that V2{q,t) = Vi{q,t). 

Let us now consider the class P3 of r^-symmetric convex polyominoes with /ii = /i2 = h. In 
this case, the added decorations on the super-hexagon can only occupy one of the sectors A, B, 
C 01 D shown in Figure 1141 with the exception of the sectors Ar\ B and C D D which can be 
simultaneously occupied. 




a) b) 

Figure 14: Fundamental regions of r^-symmetric convex polyominoes in V3 



Denote by B^ the class of admissible decorations in the sector B over an hexagon of side h, 
and by Bh{q,t) its generating series, where the variables q and t mark the added area and half 
perimeter, respectively. We have 

Bh{q, t) = E ( - (9'* + '^i,r,i{q, t)) + V2,rAQ, + E *'"'^r,i((7, t)] (47) 

r=l V j>2 / 
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The generating series will be the same for the decorations located in the sectors A, C or D, 
for symmetry reasons. However, in the term ABh{q,t), the decorations which, like the one 
shown in Figure [Tib , are located in the intersection sectors APiB or C Ci D, are counted twice. 
Observe that these decorations are in fact stack polyominoes of the class TO with generating series 
TO/j_i(t, q) — l, whence the correcting term — 2(T0ft_i(t, q) — l)- Lastly, the term {TOh~i{t, q) — l)'^ 
counts the simultaneous decorations in the sectors An B and C H D and the term 1 is added 
for the empty decoration. Globally, we obtain 

V3{q,t) = ^^3(2h-l)^3M/^-l)+l (^4Bh{q\t^)-4TOh-i{t\q^)+TOl_,{t\q^)+A) (48) 

h>l 

5.2.2 The rotation center is on a vertex 

The class of r^-symmetric convex polyominoes whose rotation center is on a vertex is denoted 
by Q. There are three cases: hi > h2, h2 > hi and hi = /12 to which correspond three series 
Qi{q^ t), Q2{q, t) and Qsiq, t) and two types of central vertices as in FigurelTKl The computations 




Figure 15: r^-symmetric convex polyominoes in Qi 



are similar to the preceding case. The decorations are placed over a pseudo-hexagon and we find 

h+l h -, 

Qiiq,t) = 2 ^ i6h,g3h,' ^ ^ ^ ^Po,r,/(9^^')^?o,.,/(g^^') + ^5o,r,/(9^^')^?l,.,/(9^^') + 

h>l 1>1 r=l s=l 

25o,r,/(9^^^)^2,.,/(g^^^) + 25l,r,/(g^^^)^?l,.,/(9^^^) + 

Vi^rAi^ ^^)Po,s,/(g^ i^) + ^^2,r,/(g^ ^^)^?o,s,/(g^ i^)) (49) 

= Q2{q,t) 

and 

Q3{q,t) = 2j2t^\^^" [A + 2Bh{q',t^)+2Bh+i{q\t')-2TOH^iit\q^) 

h>l 

-2T0h{t\ q^) + TOh-i{t\ q^)TOh{t\ q^)) . (50) 
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5.2.3 Global result 

Finally, 

\Fix{v^)y = 2Vi{q,t)+Vs{q,t) + 2Qi{q,t) + Q3{q,t). (51) 
5.3 Symmetry with respect to the tt radian rotation r"^ 

Here, the rotation center can be in the middle of an edge or of an hexagon. See Figure ITHl If the 




rotation center is in the middle of an edge, there are three similar cases corresponding to the 
three types of edges. Consider the case of the horizontal edge and denote by A, the corresponding 
class. Such a polyomino P is shown in Figure [Wk . Denote by K the left fundamental region of 
P, including the central column. Note that this column is of even length. 

If the rotation center is the middle of an hexagon, we denote by TC the corresponding class. 
In this case, the central column is of odd length. See Figure 116b . The polyominoes K which 
can occur as a fundamental region in one of these two cases are 

K = Coo + Coi + Cio + Cu + Co2 + C20. (52) 

Recall that the series Cij^ni^, q, t) is defined by equation H25|) . with the index n representing the 
size of the last column. We then have 

Aix,q,t) = 5I^fc74fcTT('^oo,2fc + 2Col,2fc + Cll,2fc + 2Co2,2fc)(2;^g^^^). (53) 

and 

n{x, q,t) = 2fc+U4fc+3 ('^00,2fc+l + 2Coi,2fc+l + Cii^2fc+l + 2Co2,2k+i){x^ , , t^)- (54) 

fc>o ^ 

Finally, 

|Fix(r3)|^_i = 3^(1, g, t) + n{l, q, t). (55) 
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6 Two generator symmetry classes 



6.1 Symmetry with respect to Dq 

Since Dq =< r, ds2 >, convex polyominoes belonging to Fix(T>Q) consist of super-hexagons with 
symmetric stack decorations (see the section l??T|) . We obtain 

|Fix(Z)6)|g,t = J2 ^'^''^-')'?"^^'^-'^+'TSO;,_l(^^ (56) 

h>l 

6.2 Symmetry with respect to F3 i =< r^, ds2 > 

We are guided by the 2tt/3 rotation symmetry class studied in section 15.21 The cases hi > /12 
and /i2 > hi are impossible because of the ds2-symmetry. There remains the case /ii = /i2 = h. 
If the rotation center is in the middle of an hexagon (case V3), the sides of the superhexagon 
are decorated by symmetric stacks of type TSO. Moreover the ds2-symmetry implies that the 
decorations are in the sectors A n B and C H D. By r^-symmetry, three of these stacks are 
identical and the three others also, whence the formula {TSOh-i{t^ , 9^))^- If the rotation center 
is a vertex (case Q3), we rather find TSOh-i{t'^,q'^)TSOh{t'^,q'^)- Consequently, 

|Fix(F3,i)|,, = ^^='(2.-l)^3/.(.-l)+l(^^0._l(^^g3))2 + 

h>l 

2 J2 t'^^^'^'TSOhMt', q')TSOh{t^ q'). (57) 

h>l 

6.3 Symmetry with respect to H31 =< r^,v > 

We refer again to section 15.21 The case where the center is a vertex is impossible because of 
the vertical symmetry. There remains the case where the center is an hexagon and the three 
subcases hi > h2, h2 > hi and hi = h2 define three subclasses TZi, 7^2 and 7^3, respectively. 
For the case where /ii > /i2, the part D* of the decoration (see Figure IT^ is in fact the mirror 
image v ■ D of D. Hence we obtain 

ni{q,t) = ^^3(2'^-l),3'^('^-l)+^5:x:4I(^o,nK^^^') + ^^l,nK'?^^^)) (58) 

h>l l>lr=l^ 

= T^2{q,t). 

If hi = h2 = h, the decorations in the sectors AOB and C f] D are mirror images of each 
other and we find 

7^3(g, t) = Y, t^'^^''~^^q^^^^~^^^^TOh_i{t^,q^). (59) 

h>l 

Finally, 

\Fix{H3,i)\g,t = 27^l((7,^) +7^3((/,^)• (60) 
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6.4 Symmetry with respect to ^2,3 =< r^, h > 

Observe that -D2,3 =< h, v >. We thus refer to sections 14.11 and 14.21 on v- and h-symmetric 
polyominoes, respectively. In order to obtain a Z)2,3-synimetric convex polyomino, it suffices 
to take a v-symmetric polyomino whose fundamental region K (see the figure lll|) is itself h- 
symmetric. 

The series CSoo(a;, <?, f, t), HSii(x, u, f , t) and CSii(a;, g, f, t) are the h-symmetric analogues 
of the series Coo(x, q, t), Hii{x,q,u,v,t) and Cu{x,q,v,t) of sections 12331 and 123 We have 

CSoo,fc(x, q, t) = b'=]CSoo(x, q, V, t) = t^'^TSkixt, q), (61) 
}lSii{x,q,u,v,t) = + xqvt"^ IlSu{x,q,u,vq,t) + 



1 — quvt^ 



xt f RSn{x, q,u,v,t) 

— E.bn[x, q,u,vq,t) - v{ )^=o 

qv \ V 



^HSn,fc(x,(?,^;,^)u^ (62) 

k>l 



CSii{x, q,v,t) = '^-^CSoo^i{x,q,t)RSii^i-i{x,q,v,t) (63) 



+2i 

i>2 ^ 

^ CSii,fc(x,g,t) 
fc>i 



and finally 



|Fix(Z?2,3) = E (cSoo,i(l, g', t') + CSii,,(l, q^,e)) . (64) 



7 Conclusion 

It is now possible to use Burnside's formula H12|) . with ^ = C, to enumerate the free convex 
polyominoes, according to area and half perimeter. Some numerical results are given in tables 
1 and 2, according to area only (up to area 20) or to half perimeter only (up to half perimeter 
16). See under the column "Orbits". 

It is also possible to enumerate asymmetric convex polyominoes with the help of formula 
H16|) . with T = C . Some results are found in the tables 1 and 2. It is clear on these tables that 
almost all convex polyominoes are asymmetric. 

All these numerical results were verified experimentally by an exhaustive computerized enu- 
meration. 
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Area 


id 


h 


V 


r 


J.2 




Orbits 


v& 


^31 


H31 


D21 


Asym 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 
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2 


3 


1 
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3 
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1 





3 


11 


3 


3 





2 


3 


3 
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4 


38 


2 


4 





2 


12 


6 








2 


2 


24 


5 


120 


6 


10 








12 


15 











2 


72 


6 
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6 


12 





6 


42 


38 





2 





2 
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7 


939 


9 


27 


1 


3 


37 


91 


1 


1 


3 


3 
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8 
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12 
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4 
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9 
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17 


63 





12 


99 
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3 
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10 
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20 


66 





6 
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2 


4 
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11 
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30 
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6 
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12 
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38 
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18 
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2 





4 
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13 
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46 
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1 


13 
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1 


1 
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14 
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62 
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69 
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V 


r 


r^ 


j,3 


Orbits 


-De 
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H31 
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1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 





4 
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1 


1 








3 


1 











1 





6 


2 


2 








2 





1 





2 











7 


12 


2 


4 








6 
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2 





8 


18 


2 














2 














12 


9 


59 


5 


9 


1 


5 


19 


11 


1 


3 


3 


3 


24 


10 
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8 














12 














96 


11 
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10 
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48 


39 











6 


204 


12 
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14 








12 





65 
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672 


13 
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25 


59 
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36 
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